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O 

It will be shown that starting from a coordinate system where the 6 
^ phase space coordinates are linearly coupled, one can go to a new coordinate 

^ system where the motion is uncoupled by means of a linear transformation. 

The original coupled coordinates and the new uncoupled coordinates are 
Q related by a 6x6 matrix, R. R will be called the decoupling matrix. It will 

O be shown that of the 36 elements of the 6x6 decoupling matrix R, only 

^ 12 elements are independent. This may be contrasted with the results for 

in motion in 4-dimensional phase space, where R has 4 independent elements. 

^ A set of equations is given from which the 12 elements of R can be computed 

^ from the one period transfer matrix. This set of equations also allows the 

^ linear parameters for the uncoupled coordinates to be computed from the 

JpH one period transfer matrix. An alternative procedure for computing the 

O linear parameters and the 12 independent elements of the decoupling matrix 

^ R is also given which depends on computing the eigenvectors of the one 

period transfer matrix. These results can be used in a tracking program, 
where the one period transfer matrix can be computed by multiplying the 
transfer matrices of all the elements in a period, to compute the linear 
parameters and the elements of the decoupling matrix R. 



* Work performed under the auspices of the U.S. Department of Energy. 
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Introduction 



1. Introduction 

It will be shown that starting from a coordinate system where the 6 phase space 
coordinates are linearly coupled, one can go to a new coordinate system, where the motion 
is uncoupled, by means of a linear transformation. The original coupled coordinates and 
the new uncoupled coordinates are related by a 6 X 6 matrix, R. R will be called the 
decoupling matrix. 

It will be shown that of the 36 elements of the 6x6 decoupling matrix i?, only 12 
elements are independent. This may be contrasted with the results-^ for motion in 4- 
dimensional phase space, where R has 4 independent elements. A set of equations is given 
from which the 12 elements of R can be computed from the one period transfer matrix. 
This set of equations also allows the linear parameters, the a^, z = 1, 3, for the uncoupled 
coordinates, to be computed from the one period transfer matrix. 

An alternative procedure for computing the linear parameters, /3i,ai, i = 1,3, and 
the 12 independent elements of the decoupling matrix R is also given which depends on 
computing the eigenvectors of the one period transfer matrix. 

These results can be used in a tracking program, where the one period transfer matrix 
can be computed by multiplying the transfer matrices of all the elements in a period, to 
compute the linear parameters ai and z = 1,3, and the elements of the decoupling 
matrix R. 

The procedure presented here for studying coupled motion in 6-dimensional phase 
space can also be applied to coupled motion in 4-dimensional phase space, where it may 
be a useful alternative procedure to the procedure presented by Edwards and Teng"^. In 
particular, it gives a simpler programing procedure for computing the beta functions and 
the emittances for coupled motion in 4-dimensional phase space. 
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2. The Decoupling Matrix, R 



The particle coordinates are assumed to be x, y, Py^ pz- The particle is acted 
upon by periodic fields that couple the 6 coordinates. The linearized equations of motion 
are assumed to be 

rl x ^ ^ 

A(s) X 



dx 
ds 



[2-1] 



X 
Px 

y 

Py 

z 

Pz 

where the 6x6 matrix A[s) is assumed to be periodic in s with the period L. Note that 
the symbol x is used to indicate both the column vector x and the first element of this 
column vector. The meaning of x should be clear from the context. This kind of double 
use of a symbol will be used in several places in this paper. The 6x6 transfer matrix 
T(s, So) obeys 



x{s) = T (s, So) X (so) 
dT 
ds 

It is assumed that the motion is symplectic so that 

TT = / 
T =ST S 



(2-2) 



(2-3) 



where / is the 6x6 identity matrix, T is the transpose of T and the 6x6 matrix S is 
given by 



f2-4a) 
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(2-4b) 
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The Decoupling Matrix, R 



where the 2x2 matrix S is given by 



S 



1 
-1 



(2-4d) 



Note that the same symbol S is used to indicate the 6x6 and 2x2 symplectic matrix. 
The meaning of S should be clear from the context. 

The 6x6 transfer matrix T(s, sq) has 36 elements. However, the number of independent 
elements is smaller because of the symplectic conditions given by Eq. (2-3). There are 
15 symplectic conditions or (A;^ — k)/2 where k = 6. This can be seen by noting that the 
equation TT = / can be written as 



(2-5) 



L> = T ST -S = 
D = -D 

Since D is anti-symmetric, it has (36 — 6)/2 = 15 independent elements. Eq. (2-5) 
then represents 15 symplectic conditions. The transfer matrix T then has 21 independent 
elements. 

One can also introduce the one period transfer matrix T(s) defined by 



T{s + L,s) 



(2-6) 



T(s) is also symplectic and has 21 independent elements. 

The 15 symplectic conditions for T can be written down more explicitly. T can be 



written in terms of the 2x2 matrices 




= 1, 


3,J = 






T12 


Ti3 


T = 


T21 


T22 


T23 




T31 


T32 


T33 



f2-7a) 



Using T =ST S, one can show that T is given by 



Tn 


T21 


T31 


T12 


T22 


T32 


Ti3 


T23 


T33 



(2-7b) 



Then TT = / gives the conditions 

^ ^ Tj^Tj^ = Sij 



I = 1,3 J = 1,3 



f2-8a) 



it=l,3 
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which can be written as 



|Tii| + IT12I + iTisI = 1 
IT21I + IT22I + IT23I = 1 
IT31I + IT32I + iTssI = 1 



(2-8b) 



T11T21 + T12T22 + T13T23 = 
T11T31 + T12T32 + T13T33 = 

T21T31 + T22T32 + T23T33 = 

To get the first 3 equations, one uses the resuh for 2x2 matrices 

where is the determinant of Tij. Eqs. (2.8) are the 15 symplectic conditions for T. 

One now goes to a new coordinate system where the particle motion is not coupled. 
The coordinates in the uncoupled coordinate system will be labeled t/;,^^. It 

is assumed that the original coupled coordinate system and the new uncoupled coordinate 
system are related by a 6 X 6 matrix R{s) 



X 


= R u 






u 






Pu 






V 


u 










Pv 






w 






-Pw. 



(2-9) 



R{s) will be called the decoupling matrix. 

One can introduce a 6 X 6 transfer matrix for the uncoupled coordinates called P(s, so) 
such that 

u{s) = P {s,so)u (2-lOa) 



and one finds that 



P(s,sn) = R-Us)T(s,so)R(so, 



(2-lOb) 



one can also introduce the one period transfer matrix P{s) defined by 

P{s) = P{s + L,s) 

P{s) = R-^{s + L)T{s)R{s) 



(2-11] 
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The Decoupling Matrix, R 



The decoupling matrix is defined as the 6x6 matrix that diagonahze P{s)^ which 
means here that when the 6x6 matrix P is written in terms of 2 X 2 matrices it has the 
form 



Pn 













P22 





(2 








As 





2-121 



P 



where Pij are 2x2 matrices. P wiU be caUed a diagonal matrix in the sense of Eq. (2-12). 

The definition given so far of the decoupling matrix i?, will be seen to not uniquely 
define R and one can add the two conditions on R that it is a symplectic matrix and it is 
a periodic matrix. Thus 

RR = I 

R=SRS (2-13) 
R{s + L) = R{s) 

The justification for the above is given by the solution found for R{s) below in section 4. 

Because T(s,so) and R{s) are symplectic, it follows that P(s,so) and P{s) are sym- 
plectic. Eq. (2-8) and (2-9) can be rewritten as 



P{s,so) = R{s) T{s,so) R{so) 
P{s) = R{s) T{s) R{s) 

It also follows that the 2x2 matrices has 3 independent elements as |Pii | 
1. Eq. (2- 14b) can be written as 



(2-14a) 
(2-14b) 



IP22I = IP33I 



T(s] 



R(s) P(s) R(s] 



(2-15) 



Eq. (2-15) shows that R has 12 independent elements, as T has 21 independent elements 
and P has 9 independent elements. As R has only 12 independent elements, one can 
suggest that R has the form, when written in terms of 2 X 2 matrices. 



R 



qil R12 Ru 

R21 q2l R23 
R31 R32 q^I 



(2-16) 



where (j'l, (j'2, qs are scalar quantities, the Rij are 2x2 matrices and / is the 2x2 identity 
matrix. The matrix in Eq. (2-16) appears to have 27 independent elements. However, R 
is symplectic and has to obey the 15 symplectic conditions, and this reduces the number 
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of independent elements to 12. The justification for assuming this form of i?, given by Eq. 
(2-16), wiU be provided below where a solution for R will be found assuming this form for 
R. 

The 15 symplectic conditions for R may be written, using Eq. (2-8b), as 

91 + \Rn\ + |i?i3| = 1 
\R2l\+ql + |i?23| = 1 
|i?3l| + |i?32| +g| = 1 

(2-17) 

qiR2i + Ri2q2 + RuR23 = 
qiRsi + RnR32 + Ruqs = 

-R2I-R3I + 52-^32 + -R23?3 = 

Using Eq. (2-16) for R and the symplectic conditions Eq. (2-17), one can, in principle, 
solve Eq. (2-15) for R and P in terms of the one period matrix T. This was done by 
Edwards and Teng"^ for motion in 4-dimensional phase space where T has 10 independent 
elements, R has 4 independent elements and P has 6 independent elements. An analytical 
solution of Eq. (2-15) for the 6-dimensional case was not found. However, a different 
procedure for finding P and R will be given below in section 4 by finding the eigenvectors 
of P, using the eigenvectors of the one period matrix, T. 

The 2x2 matrices Pn. P22, -P33 which make up P each have 3 independent elements 
and can be written in the form 



Ai = 



cos -|- «! sin /3isin?/'i 
— 1/71 sin cos — «! sin 



71 = (l + «2)//3i 



f2-18) 



with similar expressions for P22 and P33 Eq. (2-18) and the similar expressions for P22, P33 
can be used to define the three beta functions (32 and (S^. This definition of the beta 
functions will be further justified below where one finds expressions for the three emittance 
invariants ei, €2 and 63 in terms of ai, /32, 02 and /33, a^. If one could solve Eq. (2-15) 
for R and P in terms of the one period transfer matrix T, then Eq. (2-18) could be used 
to find the beta functions (3i from the one period transfer matrix. A different procedure 
for finding the linear parameters ai, /32, 02 and /33, will be given below. 
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The Linear Parameters j3 , a, and ip and the Eigenvectors of the Transfer Matrix 



3. The Linear Parameters (3, a, and tjj and the Eigenvectors of the Transfer Ma- 
trix 

In this section, the eigenvectors of the one period transfer matrix, P, wiU be found and 

expressed in terms of the hnear periodic parameters /3, a and These wiU be used below 

to compute the linear parameters from the one period transfer matrix T. They will also 

be used to find the three emittance invariants ei, 62 and 63 and to express them in terms 

of the linear parameters a^, z = 1,3. 

The uncoupled transfer matrix obeys 

d 



— = P{s,so)=B{s) P{s,so] 
as 

B = R A R+^R 

as 



(3-1] 



This follows from Eq. (2-2) and Eq. (2-14). 

One sees from Eq. (3-1) that B[s) is a periodic matrix, B[s -\- L) = B[s). It can also 
be shown that i? is a diagonal matrix in the sense that written in terms of 2 X 2 matrices, 
it has the form 



B 



Bn 
B22 
533 



f3-2l 



f3-3) 



To establish Eq. (3-2), one needs the equation for dP / ds 

dp 

— = B P - P B 

ds 

Eq. (3-3) follows from the result 

P{s)=P{s,so)P{so)P{so,s) (3-4) 

which only requires that P[s -\- sq -\- L) = P(s,so), which follows from P(s,so) = 
R{s)T[s^ so)R[so) and T[s -\- So -\- L) = T(s,so)- Eq. (3-4) when differentiated gives Eq. 
(3-3), using the results 



^P(30,3) 



B(s)P(s^Sr 



(3-5) 



-P{so,s)B{s) 

The last equation follows from P(s, so) Pi-so^ -s) = One can show that B12 = by using 
Eq. (3-3) 

Ap^2 = = P12P22 - PiiBu = B12 {P22 - Ai 

5i2 =0 



(3-6) 
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In the same way one can show that Bij = if i ^ j . It has now been estabhshed that the 
uncoupled coordinates indicated by the column vector u obeys the equations of motion 



du 
ds 



B 



which can be written as 



d 
ds 



u 




Bn 










u 


V 














V 


w 










-B33 




w 



u 



u 




V 




w 




,v = 




, w = 




_Pu_ 




_Pv _ 




_Pw _ 



f3-7a) 



(3-7b) 



The three sets of coordinates, u^v^w are uncoupled, and the u coordinates obeys 



—u = Biiu 

ds 



f3-8l 



As the 2x2 matrix Bn is periodic, one can show^ that the eigenvector of the transfer 
matrix for u is 



Ul 



.1/2 



Pi 



1/2 



exp {itpi 



f3-9a) 



_f3{'^{-ai + i} 
u-^ S Ul = 2i 

with the eigenvalue Ai = exp(z/ii). /3i(5), ai[s) are periodic functions and the phase 
function rpi = /iis/L + gi{s) where gi{s) is periodic. One may notice that the Bn matrix 
does not have the properties usually assumed in the large accelerator approximation, that 
(-Sii)ii = (-821)22 = and (-Bii)i2 = 1- Thus obtaining Eq. (3-9) requires a study^ of 
the more general equation (3-8). In particular the relationship between ai, rpi are now 
given by^ 



ds 

V'l = / (i?ll)i2 







/5: 



ai 



(3-9b) 



One can write down similar expressions for the eigenfunctions of the transfer matrix of 
V and t/J, in terms of /32, 02, V'2 and /33, a3, ?/'3, with the eigenvalues A2 = exp(z/i2), 
A3 = exp(z/i3). 
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Computing the Linear Parameters j3 , a, ip from the Transfer Matrix 



One can now write down the eigenvectors of the P matrix using Eq. (3-9). These 
eigenvectors wiU be caUed ui, U2, U3, U4, U5, ue, each of which is a 6 X 1 column vector. 



Ul 



■1/2 



1/2 



-ai + i] 



Us 












si" 

,33-'" (-03 + !) 



exp {itpi 



Us 



(3. 






1/2 



^^"^ (-02 + i) 





exp {i'tp2] 



(3-10) 



exp [itps] 



U2 = t^l, U/l = U3, Uq = 

no* no* no* 

Ul Sui = u^Suz = u^Sus = 2i 
with the eigenvalues Ai = exp(z/ii), A3 = exp(z/i2), A5 = exp(z/i3), A2 = A^, A4 = A3 and 
Ae = A5. 

These results for the eigenvectors will be used below to express the 3 emittance invari- 
ants ei, 62, £3, in terms of the linear parameters a^, rpi^ z = 1, 3, and to find a procedure 
for finding the linear parameters fii^a^, ^pi from the one period transfer matrix T. 



4. Computing the Linear Parameters (3, a, i/j from the Transfer Matrix 

An important problem in tracking studies is how to compute the linear parameters, /3, 
a, defined in section 3, from the one period transfer matrix. The one period transfer 
matrix can be found by multiplying the transfer matrices of each of the elements in a period. 
A procedure is given below for computing the linear parameters, which also computes the 
decoupling matrix R from the one period transfer matrix. 

The first step in this procedure is to compute the eigenvectors and their corresponding 
eigenvalues for the one period transfer matrix T. This can be done using one of the standard 
routines available for finding the eigenvectors of a real matrix. T is assumed to be known. 
In this case, there are 6 eigenvectors indicated by the 6 column vectors 
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and xq. Because T is a real 6x6 matrix, X2 = x'^, 2:4 = 3^3, x^ = x'^- The corresponding 
eigenvalue for xi is Ai = exp(z/ii) and the eigenvalue for X2 is = exp(z/ii). In a similar 
way, A2, A2 are the eigenvalues for X3 and 2:4, and A3, A3 are the eigenvalues for X5 and xq. 
The eigenvectors are then normalized so that 

X S x\ ^ 2z 

X3 S x^ = 2i (4-1) 

no* _ no* _ no* 

It will then follow that X2 8x2 = —2i, x^ 8x4^ = —2i and Xg Sxq = —2i. 

Let us now examine the eigenvectors of the one period transfer matrix in the uncoupled 
coordinate system, P = RT R as given by Eq. (2-14). The eigenvectors of P will be 
denoted by ui, U2, U3, U4, U5, ue- If xi is an eigenvector of T with the eigenvalue xi then 
one can see that ui given by 

ui = R Xl 

(4-2) 

Xl = R ui 

is an eigenvector of P with the same eigenvalue Ai. ui is related to the linear parameters 
(3i, ai, V'l by Eq. (3-10) 



Ul 



1/2 



Pi 



1/2 



Xl will be written as 



Xl 



-0,1 + i) 







Xl 
Pxl 

yi 

Pyl 

zi 

.Pzl 



exp [itpi 



(4-3) 



(4-4) 



Putting this result for ui into Eq. (4-2) and using the result for R given by Eq. (2-16), 
one finds 



Xl = qif^l^"^ exp (itpi) 
— 1/2 

Pxl = qiPi {-ai + i)exp{i'ijji 
Pxl/xi = {-ai + i) IfJi 



(4-5) 
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Computing the Linear Parameters j3 , a, ip from the Transfer Matrix 



Solving Eqs. (4-5) for ai, rpi one finds 

V'l = ph (xi) 

l/f3i=Im{p^i/xi) (4-6) 
ai = -fJiRe {pxilxi) 

where Im and Re stand for the imaginary and real part, and ph indicates the phase. 

Using Eq. (4-6), one can find the linear parameters ai, and from the eigenvector 
xi of T. In a similar way, one can find /32, 02, V'2 from the eigenvector 2:3, and /^s, as, 
from the eigenvector x^. 

One may note that one can also compute qi from Eq. (4-5) 

qi = \xi\lf3\'^ (4-7) 

with similar results for q2 and qz. A procedure is given below for computing the entire R 
matrix. Also, the tunes of the three normal modes can be computed from the eigenvalues 
of the 6x6 transfer matrix T, = exp(z/ii) i = 1,2,3. The [j^ are the phase shifts for a 
period, and give the fractional part of the tunes. 

Having found the a^, from the one period transfer matrix, one can write the 
decoupled eigenvectors ui, U2, U3, U4, U5, u% as 



Ul 



■1/2 



1/2 



Us 



U2 



-1/2 












1/2 



exp(z?/'i), U3 



-1/2 






1/2 



-0,2 + 1) 








exp {iil)2] 



(-as + i) 

U4 = U3, 



(4-8) 



exp 



Up 



Ur 



Using the decoupled eigenvector u^, and the eigenvectors of T, Xi one can find a solution 
for the decoupling matrix R. 

Let the 6x6 matrix X, and the 6x6 matrix U be defined as 

X = ^^"^ [xi X2 xz x/^ xs xq] 

U = ^^"^ [ui U2 uz U4 U5 uq] 



(4-9) 
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where the Xi and Ui are 6x1 column vectors. 

Since the eigenvectors are related by Xi = Rui one has the relationship 

X = RU 

U can be written in terms of 2 X 2 matrices as 



(4-10) 



U 



Uu 



Un 
U22 
t/33 



/^y^ exp 



^''^ (-«! + exp (zV'i) /?! ^''^ (-«! 



/^y^exp {-itpi) 

i) exp [-itpi] 



(4-11: 



-2i 



,1/2 



with similar expressions for U22 and t/33. 

The inverse matrix of U can be written as 














U22 

















-0^ 



-/^y^ exp 



(4-12) 



-2zl 



-1/2 



^^"^ (— «! -|- z) exp (itpi) (3^^ exp (^?/'l) 

It can be shown"^ that both X and U are symplectic matrices. Thus = U and one can 
then find for R from Eq. (4-10) 

R = XU (4-13) 

Eq. (4-13) provides an explicit solution for R and justifies the assumptions made in 
section 2 that R is symplectic, periodic and has the form given by Eq. (2-16). Since X and 
U are both symplectic, then from Eq. (4-13) R is also symplectic. Since the eigenvector 
xi has the form exp[iiis/ L)f[s)^ where f[s) is periodic with similar results for the other 
5 eigenvectors, all the phase factors in the product XU cancel and R is periodic in S. 

Equations (4-13) and (4-6) give a way to compute the parameters a^, rpi and the 
R matrix from the one period transfer matrix T = T{s -\- L,s). It can also be used for 
coupled motion in 4-dimensional phase space, and provides an alternative procedure to the 
one given by Edwards and Teng-^ found by solving Eqs. (2-15). The procedure given in 
this paper may be preferable for use in a tracking program, as it appears more simple to 
program. 
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The Three Emittance Invariants 



An interesting relationship that can be found from Eq. (4-13) is the connection between 
R{s) and R{sq). 

R{s) = X{s)U{s) 

R{s) = T{s,So)X{so)U{so)P{s,so) (4-14) 

R{s) = T{s,sq)R{sq)P{sq,s) 

Eq. (4-14) that relates R{s) and R{sq) is similar to the result T(s) = T(s, sq) T(so) T(so, s)- 
In section 2, the decoupling matrix was defined as the 6x6 matrix that diagonalizes 
P{s), when written in terms of 2 X 2 matrices, according to Eq. (2- 14b). It can be shown 
that R{s) also diagonalizes P{s,sq). One can write P(s,so) in terms of the eigenvectors 
of P{s) as^ 

P{s,so) = U{s)U{so) (4-15) 

Using the result for the decoupled eigenvectors Ui given by Eq. (4-8), one sees that P(s, so) 
is diagonalized when written in terms of 2 X 2 matrices. 



5. The Three Emittance Invariants 

Three emittance invariants will be found for linear coupled motion in 6-dimensional 
phase space. Expressions will be found for these invariants in terms of a^. A simple and 
direct way to find the emittance invariants is to use the definition of emittance suggested 
by A. Piwinski"^ for 4-dimensional motion. This is given by 

ei = \ xi S x^ (5-1) 

x is a 6 X 1 column vector representing the coordinates x, pj;, Py, pz- xi is a 6 X 1 
column vector which is an eigenvector of the one period transfer matrix T. xi is assumed 
to be normalized so that 

Xi S xi = 2i (5-2) 

One first notes that ei given by Eq. (5-1) is an invariant since xi S x is a Lagrange 
invariant as xi and x are both solutions of the equations of motion. Eq. (4-1) then 
represents an invariant which is a quadratic form in x, px^ y, Py^ pz- This result can 
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be expressed in terms of the linear parameters /3i,ai by evaluating ei in the coordinate 
system of the uncoupled coordinates. Since the uncoupled coordinates, represented by the 
column vector u, is related to x by the symplectic matrix i?, 

ei = \ui S up (5-3) 

ui is an eigenvector of the one period matrix P, and one sees that because of Eq. (2-14), 

xi = R ui (5-4) 

one can now use the result for u, given by Eq. (3-10) and find that 



1 



One can define 62 and 63 as 



and find 



71 = (l + «i)V/5i 



£2 = I 2:3 S 
£3 = I 2:5 S 



(5-6) 



(5-7) 



£2 = I2v'^ + 2a2vpy + (32pI 
63 = -f^w"^ -\- 2azwpw + (SspIj 
Eqs. (5-1) and (5-6) can be used in a tracking progress to compute the emittances ei, 62, 
63. The tracking program can compute the one period transfer matrix T by multiplying 
the transfer matrices of each element in a period. One can find the eigenfunctions of T, 
3^1, X3, 2:4, xs, xq and normalize them so that 

X S x\ — 2z 

X3 S x^ = 2i (5-8) 

For a given set of the coordinates x, px^ y, Py^ pz, then ei, 62, £3 can be computed using 
Eqs. (5-1) and (5-6). 

One can also show that the volume in phase space enclosed by these three emittances 
is given by 

dx dpxdydpydzdpz = 7r^eie2e3 (5-9) 
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The volume in phase space can be computed by going to the coordinate system of the 
uncoupled coordinates u, p^, p^, pw and using \R\ = 1. 

6. Summary of Results for Use in a Program to Compute the Linear Parame- 
ters and the Decouphng Matrix 

The results given above can be used to compute the linear parameters /3, a and u and 
the decoupling matrix R from the one period transfer matrix T. This section will give 
the steps required to do this using the results given above for particle motion in either 
4-dimensional phase space or 6-dimensional phase space. 

It is assumed that the one period transfer matrix T is known, having been computed 
by multiplying the transfer matrices of the elements in one period. 

6.1. Computing the Linear Parameters f3, a, v 

First, compute the eigenvectors and the corresponding eigenvalues of the one period 
transfer matrix T. To do this one uses one of the standard routines available for computing 
the eigenvalues of a real matrix. One such set of routines is the routine SGEEV (single 
precision) and DGEEV (double precision) available in the IBM Engineering and Scientific 
Subroutine Library (ESSL). 

Let us denote the eigenvectors by Xi and the eigenvalues by A^. The eigenvectors and 
eigenvalues can be arranged in pairs such that X2 = x\ and A2 = A^, and 2:4 = Xg, A4 = A3, 
and = Xg, Ae = Ag. Choose x\ to have the eigenvalue exp(z/ii) and X2 to have the 
eigenvalue exp(— z/ii) and similarly for 2:3, x/^ and 

Normalize xi, 2:3, X5 so that these 3 eigenvectors satisfy 

S Xi = 2i (6.1-1) 

The linear parameters ai can be computed from the normalized xi using Eq. (4-6) 

l/(3i = Im (pxi/xi) 

(6.1-2) 

ai = -(3iRe (pxi/xi) 
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xi and pxi are the first two elements of the column vector xi. The fractional part of the 
normal mode tune, ui can be computed from the phase advance in a period, /ii. In the 
same way one can compute /32, 02, U2 from 2:3, /is and /^s, as, U3 from 2:5, 

One may note that the normal mode tunes can be computed directly, without finding 
the eigenvectors, using the results of reference 6. This may provide a good check of the 
program. 

The diagonal elements of the decoupling matrix R qi, (j'2, qz [see Eq. (2-16)] can be 
found using Eq. (4-2) 



6.2. Computing the Decoupling Matrix R 

To compute the decoupling matrix i?, from the X matrix defined by Eq. (4-9) 



qi = \xi\/(3- 



1/2 





Form the U matrix given by Eq. (4-12) 



U 



Un 





_0 

U22 






_0 




The Uii are 2x2 matrices where 



Pi 



1/2 



(-«! - i)exp {-tpi) 
2 

(— «! -|- i) exp {itpi) 



/^y^ exp {-lipi) 
/^y^ exp (zV'i) 



1/2 



t^ll 



Pi 



1/2 




with similar results for U22 and ?7ss- 



Having formed the matrices X and U one can compute R from 



R = X U 




One may note that above results can be used for particle motion in 4-dimensional phase, 
or in 6-dimensional phase space with some minor modifications. 
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6.3. Computing the Emittance Invariants 

The emittance invariants for a set of coordinates given by the column vector x can be 
computed using Eq. (5-1) 



where the column vector xi is an eigenvector of the one period transfer matrix T. Similar 
results are found for the other invariants 62, £3. For 62 use the eigenvector 2:3, and for 63 
use the eigenvector X5. 
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